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We study freeze out process of particles across a three dimensional space-time hypersurface with 
space-like normal. The problem of negative contribution is discussed with respect to conservation 
laws, and a simple and practical new one-particle distribution for the post FO side is introduced, 
the Cancelling Jiittner (CJ) distribution. 
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I. INTRODUCTION 

Relativistic heavy ion collisions are non-equilibrium 
relativistic many body systems that can be described by 
different models, often based on the kinetic theory. The 
kinetic theory establishes relationship between macro- 
scopic and microscopic matter properties mostly by us- 
ing a one-particle distribution function. Kinetic theories 
are able to describe dilute weakly interacting systems, 
thus expanding systems where the constituent particles 
gradually loose contact. Here to simplify such a process 
we will assume a freeze out (FO) hypersurface in space- 
time, which can have either space- or time-like normal 
vector. We will use conservation laws of fluid dynamics 
to conserve energy, momentum and particle flow across 
FO hypersurface. We assume thermal equilibrium in pre 
FO side where the matter is (or can be) in quark gluon 
plasma (QGP) phase. For the pre FO side the Bag Model 
equation of state will be used which allows for supercool- 
ing and rapid hadronization afterwards. For the frozen 
out matter (post FO side) newly formed hadrons are con- 
sidered, which are non interacting with each other, and 
they are described by a one-particle distribution func- 
tion. The invariant number of conserved particles (world 
lines) crossing a surface element, da^, of the FO hyper- 
surface is: dN — N^da^. Thus, the total number of all 
particles crossing the FO hypersurface is N = J s N^da^. 
From kinetic definition of the baryon current four vector, 
N 11 , we have: 



J —^P^fFoi^p; T,n,M p 



(1) 



Inserting it into the equation for the total number of 
particles leads to the Cooper-Frye formula pj: 



E ^P = I fFo(r,P;T,n,u^da, 



(2) 



where fFoir^p^T^n^u^ 1 ) is the unknown post FO phase 
space distribution of frozen out particles. The problem 



is to choose it's form correctly, and to determine its pa- 
rameters which satisfy all conservation laws! This was 
not done correctly in the original Cooper-Frye descrip- 
tion 0. The most used distribution for the time- like 
normal is the Jiittner distribution 0] (also called rela- 
tivistic Boltzmann distribution): 



^Juttner 



(r,p) = 



{2irh) : 



■ exp 



T(r) 



(3) 



'Electronic address: karolis@fi.uib.no 
^Electronic address: csernai@fi.uib.no 



For the time-like case p M and da^ in the Cooper-Frye 
formula are both time-like vectors, thus p^da^ > 0, and 
the integrand of integral ijSJ) is always positive. For the 
space-like normal vector case, p^da^ can be both posi- 
tive and negative. This situation is a problem, because 
the integrand in integral ma y change sign, and this 
indicates, that part of distribution contributes to a neg- 
ative current, going back into the front, while the other 
part is coming out of the front. The improvement of the 
distribution function @ was done by introducing the 
cut- Jiittner distribution Q, which makes p^da^ to 
be non-negative in Jiittner distribution by multiplying 
it with step function: fCut-Juttner = Q( J) ^da ll )f Juttner . 
The cut- Jiittner distribution has solved the above men- 
tioned problem formally, but the lack of a real physical 
solution persisted. The cut- Jiittner distribution has an 
unphysical form: the distribution is sharply cut off which 
is mathematically satisfactory, but it is hardly possible to 
imagine a realistic physical process producing such a dis- 
tribution. So, taking into account that this distribution 
describes physical particles - it has to be improved. A 
solution to this problem in kinetic theory was presented 
in ref. p|. 

In the present work we will present a new, simple dis- 
tribution, called the Cancelling Jiittner (CJ) distribu- 
tion. It solves the problem of negative contributions in 
the Cooper-Frye formula, and it has a smooth physically 
realistic form. We worked out simple relations, which 
can be used for the calculation of measurables and can 
be compared with experimental results. 
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FIG. 1: Schematic view of making the Cancelling 
Jiittner distribution: in the Rest Frame of the Front 
two Jiittner distributions (full curves) - one with neg- 
ative velocity, another with positive are subtracted. 
Dashed curve is the CJ distribution. 



II. THE CANCELLING JUTTNER 
DISTRIBUTION 

The Cancelling Jiittner (CJ) distribution, fcj, is de- 
fined by subtracting the ordinary Jiittner distribution @ 
with negative velocity, —v, from original Jiittner distri- 
bution, and and multiplying the obtained result with the 
step function (figure^: 



£ ( £ Juttner 

jcj = [Jr 



jjuttner\ 



exp 



T 



exp^pi , (4) 



where — (7, jv, 0, 0) and = (7, — jv, 0, 0) in the 
rest frame of the front (RFF). The velocity parameter, v, 
of the CJ distribution is restricted to be positive. When 
p^da a = the function vanishes at the front, even with- 
out step function 0. The role of the Q(p^da u ) part is 
just to eliminate the negative part of the distibution. 

We chose this construction from two Juttner distribu- 
tions because: 1) It automatically includes the cut, and 
furthermore a smooth but rapid cut, thus the unrealis- 
tic profile of the earlier proposed cut- Juttner distribution 
is not present. 2) It resembles the distribution we ob- 
tained from a kinetic model quite well, in any case much 
better than the cut- Juttner distribution. 3) The formu- 
lation is still analytic and not more complicated than the 
one arising from the cut- Juttner distribution. 

Calculations and results will be made in the rest frame 
of the front (RFF), where the surface element is da^ = 
(0, 1, 0, 0)d<x. In this frame particles cannot propagate 
across the front in the negative x direction. We assume 
that the pre FO side is in thermal equilibrium. If the 
particle has passed the freeze out layer it cannot scatter 
back. This layer is idealized as a front or hypersurface. 



FIG. 2: The post FO CJ distribution, f C j, in the 
rest frame of the front (RFF). Baryon mass, m — 1 
GeV, temperature, T = 100 MeV, density, n = 0.17 
fm -3 and velocity, v = 0.3 c in 'A', v — 0.5 c in 'B'. 
The values of each contour are written on the curves. 
The CJ distribution resembles strongly the post FO 
distribution obtained in kinetic theory, fig. 2 and 3 
in ref. 



In figure^] the schematic way of making the CJ distribu- 
tion is shown. The distribution obtained this way (figure 
El resembles strongly the one obtained numerically in a 
kinetic FO model [5(. In the following we demonstrate 
the use of the Cancelling Juttner distribution. 



III. CONSERVATION LAWS OF FLUID 
DYNAMICS 

The FO hypersurface is considered to be a discontinu- 
ity in space-time, thus instead of derivatives of the dis- 
continuous variables we consider the changes explicitely. 
We will use the conservation laws of fluid dynamics, but 
instead of the differential form of continuity equations we 
have to use comutator form: 



[N^dcTf,] = and [T^d(T u } = 0, 



(5) 



where [A] — A\ — Aq, A% is a post FO quantity, and 
Ao is the pre FO quantity. Moreover, entropy across FO 
hypersurface must not decrease: 



[S^] > 0. 



(6) 



This condition does not provide additional information 
about the freezing out matter, but it must be satisfied. 

Now our task is to find out the expressions for baryon 
current, energy- momentum tensor and entropy current in 
the post FO side. Having these expressions, we will insert 
them to equations of conservation laws (J5J, and then we 
will have a relationship between matter properties in pre 
and post FO sides. This calculation is straightforward 
for time-like FO where both the pre- and post FO sides 
can be characterized as perfect fluids with well defined 



3 



equation of state and energy- momentum tensor (see 
In the case of space-like FO with out of equilibrated post 
FO distribution this calculation has been demonstrated 
so far only for the cut-Jiittner distribution 

The post FO baryon current, N v , can be calculated 
inserting the post FO distribution J3J into equation QJ. 
The energy- momentum tensor, T MI/ , can be obtained in 
the same way, i.e., inserting the post FO distribution Q 
in to the definition: 



d 3 P 



f^/Fo(f,P;T,n,^) 



(7) 



However, we do not have to calculate integral QJ nor 
0, but instead we take the already calculated N^, 
and for the post FO Cut-Jiittner distribution in the 
reference frame of the gas (RFG), from Q: 
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Nqj, T^ V j and S CJ for finite mass m, in RFF are: 
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a 2 ^ 2 (a)w(l - + a 2 ii:i(a)t; 



where A, = iQ(a) - /C, (a, b), i = 1,2 1 . All other 
components vanish. From these expressions, one obtains 
properties of the matter on the post FO side of the FO 
front. 



IV. CALCULATING THE RELATIONSHIP OF 
MATTER PROPERTIES 

Now we will study in more detail the change of mat- 
ter properties across the FO hypersurface with space- 
like normal, having local equilibrium in pre FO side, and 
non-equilibrated matter on the post FO side. In order 
to connect both sides of the FO hypersurface we will use 
conservation laws of hydrodynamics. 

We are taking baryon current and energy momentum 
tensor for the cut-Jiittner distribution in the RFG with 
mass zero limit 2 . Such a approximation is possible when 
energies are high compared to the masses of particles. 
From Q Nft FG and T^ FG read as: 



AT 



1 



where j = sign(u), h — TrT 3 e^/ T {ith) 3 , [i is chem- 
ical potential of the pre FO matter, a = m/T, b = 
a/y/l - v 2 , v = da /da x , A = (2 + 2b + b 2 )e~ b , 
(l + b + b 2 /2 + b 3 /6)e- b and 



B 



2 n (nl) 
fcn(a,b) = a 7 

(2n)\ 



dx{. 
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To get the post FO baryon current and energy- 
momentum tensor for the CJ distribution in the rest 



i.e. K n (a) = K n (a,a). All other components of T 7 '", N 11 
and 5* M vanish. 

The derivation of N GJ and T^j will be presented in 
detail in the mass zero limit in the next section. However, 
the way of deriving these quantities is the same for both 
cases: with mass zero limit, or without. The results of 



0, all components of T^" and JV equal 



1 when v = 0, i.e. 
zero. 

2 notations of baryon current, energy-momentum tensor and en- 
tropy current with mass zero limit and without are the same; 
from now on m = is assumed. 
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frame of the front (RFF), we are making a Lorentz trans- 
formation from the RFG to the RFF of the components 
of the cut-Jiittner distribution. Thus, we get Nj^ FF and 
Topp for the cut-Jiittner distribution in mass zero limit: 

N£p F = ^-(-v 3 + 3v + 2,v 2 + 2v + 1,0,0) 



give an identical result, because the CJ distribution is not 
an equilibrium distribution. In further calculation let us 
use Eckart's definition of the flow. Having expression 
for the flow, we can obtain other macroscopic quantities. 
Invariant scalar density n is: 

n = N£ju*. (13) 



Energy density, e, and entropy density, s, reads as: 
?Tff = l^(v 5 -6v 3 +2v 2 + l2v+6), e s = u E T^u E , s E = S»u E . 
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RFG 



The energy-momentum tensor components T yv and T zz 
are not changing, because our reference frame is chosen in 
such a way, that the Lorentz transformation is influencing 
only time and one spatial component x. All other energy- 
momentum tensor components are equal to zero. 

Any of the components presented above for the CJ 
distribution, for example, Nqj in the RFF frame, can be 
found in such a way: 

NcM = Kff{v) - N° RFF (-v). 

This follows from the definition of the CJ distribution 
function, eq. Thus, the components of the baryon 

current of the CJ distribution take the form: 
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FIG. 3: The Eckart's flow velocity, it, of the frozen- 
out particles, described by CJ distribution in RFF, 
versus velocity parameter of the distribution, v. 



3v 



,v,0,0 



(8) 



and the components of the energy-momentum tensor: 
T°° = »T7V-3v 3 +6«), — T^°j — 4hT^ 2 v, (9) 



TS x J = nT 1 2 (v 3 +3v), T™=T^j = nT 



-v d + 3v 



(10) 

Now we can introduce flow, particle density, energy den- 
sity and entropy density for the frozen-out matter. The 
four- flow of the particles, u^, can be described in different 
ways: i) Eckart's definition: 



u 



flow.E 



iV CJ 



y/NZjN VtC j 



(11) 



where the flow u E is tied to conserved particles, ii) Lan- 
dau's definition: 



L flow,L 



V U P T CJ U " 



(12) 



where the flow is tied to energy flow. In case of Can- 
celling Juttner distribution there two definitions do not 




v[1/c] 



FIG. 4: The proper density of the frozen-out parti- 
cles, n, described by CJ distribution, versus the veloc- 
ity parameter of the distribution, v, with fixed chem- 
ical potential, /i = lQMeV , and temperature param- 
eter, T = VOQMeV. 

As we can see in figure |3 where the flow velocity is 
plotted versus velocity parameter of the C J distribution, 
the Eckart's flow has the minimum value 0.66c in the 
mass zero limit. Figure 0] shows the post FO density 
,n E , dependece on velocity parameter, v. 
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V. FREEZE-OUT FROM QGP 

We must find expressions for the N{f and T^ v in the pre 
FO side to use conservation laws (J3J) across the surface 
element da^. For perfect fluids the energy-momentum 
tensor (in any reference frame) can be written as 0: 

T" v = (e + p)u"u v -pg"", (14) 



The Bag Model Equation of State (EoS) for quark 
gluon plasma is used to describe the pre FO matter state. 
It is assumed that quarks and gluons exist in perturba- 
tive vacuum, where plasma contains 2(iV ( ? — 1) gluons 
and 2N c Nf quarks (N c and Nf are the number of colors 
and number of flavors) @ . The Bag Model EoS is based 
on Stefan-Boltzmann EoS, including a bag constant B. 
From [6( pre FO quantities are expressed as: 



where e is the energy density and p is pressure in the local 
rest (LR) frame. For the baryon current, we use Eckart's 



definition for flow: 
definition follows: 



N (o) =n o u (oy 



From this 



(15) 



eo 



n 



1 



54tt 2 ' 



A 4 b 



(he) 3 



l 



(he) 



where u^y. = (1,0,0,0) and the invariant scalar, uq — 

J N^N^ , is the particle density. 

To determine the parameters of the CJ distribution 
we will use the conservation laws introduced by eq. 
(JHJ. By using the conservation law for the four-current 



[N» d& f_i\ = N^jda^ - N^da^ = 0, where daf' = 
7o(^o, 1,0,0) in the pre FO local rest frame, and da^ = 
(0,1,0,0) in the RFF, we get the equation: 



hjv - n Q j v = 0. 



(16) 



Note that post and pre FO expressions have to be in the 
same reference frame, in order to use conservation laws 
correctly. Nevertheless, N{f and da^ do not have to 
be transformed to the RFF, because the baryon number 
crossing the surface is an invariant scalar. However, the 
energy-momentum current, T^da^, is not an invariant 
scalar, thus it must be Lorentz transformed to the RFF. 
After the transformation T^ u gets a form: 



(0)RFF~ 



7o( e o+"oPo) la(-eov ~pov Q ) 

7o(-e Q vo-p v ) 7o( e o«o+Po) 

po 

po 



Using the conservation laws for the energy-momentum 
tensor, [T M " da^] = T^jda^ — T^dcr^ = 0, we obtain 
two more equations: 



47 nTv - 7 (v e + v po) = 



(17) 



7 2 nT(i; 3 + 3v) - 7o 2 (ufa + Po ) = 0. (18) 

Entropy should not decrease, [S^dcr^] — S^jda^ 
S^da^ > 0, it leads to the condition: 



n 7 w (l - ip) ~ so7o«o > 0, 
which must be satisfied. 



(19) 



Po 



where B 



A 



(he) 



Here /x is baryon chemical potential. 

Using above and equations (|16l - I18fl we can evaluate 
how post FO matter properties depend on the pre FO 
side matter properties. It is enough to fix four quantities 
in the Bag Model EoS to obtain the properties of the 
post FO matter. To describe the pre FO side we will 
use: initial temperature, To, bag constant, As, initial 
baryon density, no, and initial velocity, vq, to have pre 
FO energy density, pressure, entropy density and baryon 
chemical potential. 

Nevertheless, the CJ distribution has a negative as- 
pect: not for all initial values it is possible to calcu- 
late post FO matter parameters. There are two rea- 
sons: 1. Entropy must not decrease, 2. The maxi- 
mum of the Juttner distribution function must be on the 
positive velocity side. When maximum of the Juttner 
distribution is at velocity zero (fj u ttner(v — 0)\rff = 
max fjuttner(v)\RFF), the CJ distribution is a 'zero' func- 
tion (fcF(v)=0,y v). This problem depends on how pre 
FO properties such as velocity, vq, density, no, bag con- 
stant, As, and initial temperature, To, are chosen. 

The boundary conditions for pre FO velocity, Vq and 
density, no, with fixed bag constant, As, and initial tem- 
perature, To j are shown in figure Such a conditions 
for the pre FO matter, which are to the right from the 
curves in figure[3 have to be satisfied for the CJ distribu- 
tion because we are dealing with a FO hypersurface with 
space-like normal. The shape of the curves are influenced 
by the way of making the CJ distribution i.e. when the 
post FO velocity becomes imaginary. The cutoff of the 
curves is influenced by the entropy condition. In order to 
have perspicuity we compare two curves: one with initial 
temperature T = 60 MeV, another with T = 100 MeV 
(Ab = 200 MeV for both curves) . Having higher temper- 
ature, the density (or velocity) can be lower than in the 
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FIG. 5: The boundaries of the pre FO density, no, 
and velocity, vo , to have correct post FO distribution 
function, fcj, in the rest frame of the front (RFF) 
with fixed initial temperature, To, and bag constant, 
As. The CJ distribution is applicable for initial con- 
ditions which are to the right from the curves, i.e. for 
large flow velocities in RFF. 



FIG. 6: Change of the velocity parameter across the 
FO hypersurface, in RFF, described by the CJ distri- 
bution for the different initial parameters: 1. no = 0.5 
mi -3 , A = 150 MeV full curve; 2. n = 1 mi -3 , 
A = 200 MeV dashed curve; 3. n = 0.2 fm" 3 , 
A = 150 MeV dotted curve and T = 60 MeV for all 
curves. 



lower temperature case. On the other hand, to insure the 
entropy condition, the density in the higher temperature 
case must be higher (more than 0.5 fm -3 ) than in the 
case of lower temperature, where it must be just more 
than 0.1 fm -3 . Temperature is influencing the length of 
the curve and moves it to the left when temperature in- 
creases, and bag constant is changing the gradient of the 
curve in the vo,tlq plane. 

Now, having boundary conditions we can start cal- 
culating matter properties of the post FO side. From 
equations I|16l - I18fl we get one equation for the post FO 
velocity parameter: 

v = ,/ 4( f 0+ ^-3. 
V v (eo+po) 

Note, that post FO velocity, v, temperature, T, and den- 
sity, n, are not physical quantities - they are parameters 
of the CJ distribution. As we see in figure the post FO 
velocity parameter is all the time smaller than the pre 
FO velocity. 

The post FO flow velocity, Ufi ow , is calculated using 
Eckart's definition of the flow using formula ljll|l . Figure 
shows the post FO baryon flow velocity dependence on 
the pre FO velocity, v - We observe, that the flow does 
not decrease in the same way as the post FO velocity 
parameter, v. Only in the case of low initial density 
(rio = 0.2 fm" 3 in figured) we have decrease of velocity 
going from pre to post FO side. To calculate the final 
baryon density we are using the equation ljl3|l . which 
can be rewritten for the case of Eckart's flow: 

n=^N£ J N fl<C j. 
From the results presented in the figure |S] it is seen that 
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FIG. 7: The dependence of the post FO baryon flow 
on the pre FO velocity for the different initial param- 
eters: 1. n = 0.5 fm" 3 , A = 150 MeV full curve; 
2. no = 1 fm" 3 , A = 200 MeV dashed curve; 3. 
n = 0.2 fm" 3 , A = 150 MeV dotted curve and 
To = 60 MeV for all curves. 

the baryon charge density in the post FO side is decreas- 
ing compared to the density in the pre FO side. The two 
lines are parallel to each other if initial velocity, vq, and 
initial temperature, To, are the same. 



VI. CONCLUSIONS 

We have studied the problem of negative contribution 
in the Cooper-Frye @ formula for the 3-dimentional 
hypersurface with space-like normal. The Cancelling 
Juttner distribution function was suggested as a solu- 
tion for this problem. We have showed the applicability 
and properties of the CJ distribution by using the Bag 
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0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

FIG. 8: The change of the baryon density across the 
FO hypersurface for different initial parameters: 1. 
v = 0.6 c, A = 150 MeV full curve; 2. v = 0.6 c, 
A = 200 MeV dashed curve; 3. v = 0.5 c, A = 200 
MeV dotted curve and To = 60 MeV for all curves. 



Model equation of state for QGP and conservation laws 
of hydrodynamics across the FO hypersurface. 

The C J distribution solves the problem of negative con- 
tribution in the Cooper-Frye formula and has a physical 
form unlike the cut-Juttner distribution. In future stud- 
ies one should analyze, if the conditions of applicability of 
CJ distribution are satisfied in space-like hypersurfaces 
obtained from hydrodynamical calculations. 

It will be interesting to see how large corrections this 
procedure will give to the results obtained from simple 
Cooper-Frye description. 
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